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ABSTRACT 



Context. The luminous material in clusters of galaxies exists in two forms: the visible galaxies and the X-ray emitting 
intra-cluster medium. The hot intra-cluster gas is the major observed baryonic component of clusters, about six times more 
massive than the stellar component. The mass contained within visible galaxies is approximately 3% of the dynamical mass. 
Aims. Our aim was to analyze both baryonic components, combining X-ray and optical data of a sample of five galaxy clusters 
(Abell 496, 1689, 2050, 2631 and 2667), within the redshift range 0.03 < z < 0.3. We determined the contribution of stars in 
galaxies and the intracluster medium to the total baryon budget. 

Methods. We used public XMM-Newton data to determine the gas mass and to obtain the X-ray substructures. Using the 
optical counterparts from SDSS or CFHT we determined the stellar contribution. 

Results. We examine the relative contribution of galaxies, intra-cluster light and intra-cluster medium to baryon budget in 
clusters through the stellar- to- gas mass ratio, estimated with recent data. We find that the stellar-to-gas mass ratio within 
rsoo (the radius within which the mean cluster density exceeds the critical density by a factor of 500), is anti-correlated with 
the ICM temperature, which range from 24% to 6% while the temperature ranges from 4.0 to 8.3 keV. This indicates that less 
massive cold clusters are more prolific star forming environments than massive hot clusters. 

Key words, galaxies: clusters: general- X-ray: galaxies: cluster-galaxies: luminosity function, mass function 
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1. Introduction 

Galaxy clusters occupy an unique position in the hierarchi- 
cal scenario of structure formation as they are the largest 
bound and relaxed structures to form in the Universe (e.g., 
IXozzil 2QQ7t ). An important step in understanding galaxy 
clusters is to t a ke int o account all of their components. 
Fukugita et al.l ([l998l) presented an estimate of the total 



budget of baryons in all states: most of the baryons are 
still in the form of ionized gas and the stars and rem- 
nants represe nt only 17% of the baryons. In a similar work 
Ettori ( 2QQ3I ) presented the cluster "baryonic pie" of which 



Lin et al. alyzed a sample of 13 clusters, find- 

ing that the gas-to-stellar mass ratio increases from 5.9 
to 10.4 from low- to high-mass clusters and their best-fit 
correlation for M500 — ^500 differs from L500 oc M500 by 
3<j, indicating a mass-to-light ratio increasing with mass. 
These two results strongly suggest a decrease of star for- 
mation efficiency in m ore massive environments. Recently, 



Gonzalez et al.l (|2QQ7l ) , who also used gas masses from the 



literature, estimated the contribution of stars in galaxies, 
intracluster stars and gas to the baryon budget. These au- 
thors, who took into account the intracluster light (ICL) 
to the baryon budget, confirming the previous trend of in- 
cre asing MstarZ-^gas w it h decreasi n g tem p erature found 



(|2003f ). iRoussel et al 



approximately 70% is composed of the hot intra-cluster 
medium and almost 13% composed of stars. 

The ratio of gas to stellar mass is particularly interest- 
ing because it is related to the star for mation effici ency of 
clusters. Studying the Hydra cluster, iDavid et al.l ([l99Q) 
computed the ratio of gas to stellar mass and it was found 
that the gas mass is almost 4 times that of stars. Using 
Einstein results from the literature, they extended this 

analysis to systems ranging from poor groups to rich clus- , _ ^ 

ters and a correlation between Mgas/Mgtars and the gas 'Ettoril 2QQ3I : Gonzalez et al.l 2QQ7I ) there are good reasons 
temperature was suggested. to revisit it. The most important reason is the stellar con- 



by 'D avid et all (|l99Qh : lUin et al 
(2000), found that a typical stellar contribution to the 
baryonic mass is between 5% and 20%, inside the virial ra- 
dius, even though they claim that the stellar-to-gas mass 
ratio is roughly independent of temperature. 

Although t he baryon budget h as been addressed in - 



many studies (David et all 199QI : Fukugita et al. 1998 : 
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tent contribution to the baryon fraction. Second, cosmo- 
logical simulations can precisely predict the baryon frac- 
tion within rsoo (see Sect. [7]), which is approximately the 
radius in which the X-ray observations are reliable deter- 
mined. 

To readdress this question, we choose five XMM- 
Newton public data clusters with SDSS or CFHT 
{Canada- France- Hawaii Telescope) counterparts to study 
the gaseous and the stellar baryon budgets. 

X-ray studies of galaxy clusters are particularly rele- 
vant to this issue as they provide the determination of the 
gas mass. An important task is to describe well the surface 
brightness profile, since the gas mass determination relies 
on it. 

The baryonic stellar component can be estimated by 
integrating the luminosity function to obtain the total lu- 
minosity and then adding the intra-cluster light (ICL) con- 
tribution. Using an appropriate stellar mass-to-light ratio 
we have converted the total luminosity into stellar masses. 

The paper is organized as follows: Observations and 
data treatment are described in Sect. [2l we present the 
gas density distribution in Sect. [3l the temperature spec- 
tral analysis appears in Sect. [H the dynamical analysis 
based on temperature and substructure maps is presented 
in Sect. [S) the stellar baryonic determination is presented 
in Sect. [6) the mass determination is described in Sect. [71 
Finally we present our discussion in Sect. [8] and we con- 
clude in Sect. [9l 

All distance-dependent quantities are derived assum- 
ing the Hubble constant Hq = 70 km s~^ Mpc~^, Qm = 
0.3 and, ftA= 0.7. Magnitudes are given in the AB system 
unless otherwise stated. All errors are relative to a 68% 
confidence level. 



2. Observations and data treatment 

2.1. Sample 

The objects in our sample were drawn from a set of Abell 
clusters with both XMM- Newton public archive and SDSS 
optical counterparts. We have also chosen only clusters 
with redshifts in the range 0.03 < z < 0.3 so that the 
cluster image could be well resolved and were smaller than 
the detector field of view. We thus had 6 clusters in this 
selection. Then we imposed that clusters have exposure 
times of at least 10 ks after data reduction and filtering 
(pipeline), and flare cleaning (see Sect. 12. 2p . We ended 
up with 3 clusters: 1689, A2050 and A2631. We have not 
used either the temperature or the object morphology as 
a criterion for cluster selection. 

To enlarge this sample, we added the cluster A496 that 
was observed with the Canada- France- Hawaii Telescope 
(CFHT) a nd for A2667 we used total luminosity from the 
literature (jCovo ne et al."2006bV Properties of the whole 
sample are presented in Table [H 



2.2. X-rays 

We used data from ah EPIC cameras (MOSl, M0S2 
and PN). For data reduction we used the XMM- Newton 
Science Analysis System (SAS) v6.0.4 and calibration 
database with all updates available prior to February, 
2006. The initial data screening was applied using recom- 
mended sets of event patterns, 0-12 and 0-4 for the MOS 
and PN cameras, respectively. 

The light curves are not constant and large variations 
in intensity are visible. These variations, which we call 
flares, are caused by soft energy protons produced by solar 
activity. It is better to discard flare periods and reduce 
the effective exposure time to improve the signal-to-noise 
ratio. The light curves in the energy range of [1-10] keV 
were filtered to reject periods of high background. The 
cleaned light curves exhibited stable mean count rates and 
exposure times are given in Tabled We considered events 
inside the field-of-view (FOV) and we excluded all bad 
pixels. 

In periods free of background flaring, the background 
is dominated by X-rays at low energies and particles at 
high energies. The X-ray component includes a signif- 
icant contribution from Galactic emission which varies 
with position on the sky. Since cluster emission decreases 
with distance from the center, the background compo- 
nent becomes more important towards the outskirts. The 
background was taken into account by extracting MOSl, 
M0S2 and PN spectra from the publicly available EPIC 



blank sky templates described bv Tumb et al.l (|2002 i). 



The background was normalized using a spectrum ob- 
tained in an annulus (between 9-11 arcmin) where the 
cluster emission is no longer detected. A normalized spec- 
trum was then subtracted, yielding a residual spectrum as 
shown in Fig. [TJ 

As A496 is a nearby cluster its emissivity extends up to 
the outskirts of the image and for this reason we could not 
apply this technique of background normalization. For this 
cluster, we subtracted the background without normaliz- 
ing it. If the background had not been well subtracted we 
would have had an emission excess at the outskirts. Due 
to this excess, we would have obtained from the fit of the 
temperature profile higher temperatures toward the out- 
skirts. Even so, our temper ature profile is in ag reement 
with the results obtained bv lTamura et al.l (|200ll ). 



2.3. Optical: SDSS and CFHT 

In order to investigate the optical component of A1689, 
A2050 and A2631, we downloa ded from the Sloan Dig;ital 
Sky Survey Data release 5 (Adelm an-McCarthv et al 



20071 ) the model magnitudes dereddened for all galaxies 
inside rsoo (derived in this work from an X-ray analysis, 
see Sect.[7|). Since extinction is less in the i' band we used 
this band to construct luminosity functions. The galaxy 
catalog is essentially complete down to 23.5 i' magnitude. 

A496 was observed with the Canada-France-Hawaii 
Telescope (CFHT) with the Megacam camera in 2003 
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Table 1. General cluster properties. Column (1): cluster name; Col. (2): rig ht ascension; Col. (3):dec lination; Col. (4): 
redshift; Col. (5): luminosity distance; Col. (6): hydrogen column density ( Dickev fc Lockmanl 199Qh : Col. (7): X-ray 
emission mean temperature derived in this work; see Sect. 14. 2( Col. (7): rsoo values were derived in this work. 



Cluster 


R.A 


DEC 


z 




nH 


<kT> 


rsoo 




(J2000) 


(J2000) 




(/ifo' Mpc) 


(102°cm-2) 


(keV) 


(/ifo^ kpc) 


A496 


04 33 37.1 


-13 14 46 


0.033 


144.902 


4.45 


3.96 ± 0.023 


1480 


A1689 


13 11 34.2 


-01 21 56 


0.1823 


888.870 


1.84 


8.34 ± 0.32 


2172 


A2050 


15 16 21.6 


+00 05 59 


0.1183 


551.167 


4.67 


5.85 ± 0.11 


1785 


A2631 


23 37 39.7 


+00 17 37 


0.273 


1393.23 


3.82 


6.89 ± 0.34 


2153 


A2667 


23 51 47.1 


-26 00 18 


0.23 


1146.59 


1.64 


6.56 ± 0.13 


1976 




- t * "+<4^*y I ,M 



2 ' ^ 
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Fig. 1. Comparison between observation background and Lumb background for MOSl (left panel), M0S2 (middle 
panel) and PN (right panel). The black points correspond to the observation data, the red points to Lumb scaled 
background and the green points represent the residual. 



(program 03BF12, P.I. V. Cayatte). Megacam covers a 
field of 1 X 1 square degrees with a pixel size of 0.187 
arcsec. Images were obtained in the li*, g' ^ r' and i' . 
The Terapix data center at lAP (France) reduced the im- 
ages, correcting bias and flat field and calibrating pho- 



tometrically. For more details about data reduction see 



Boue et al.l (|2QQ8[ ). 
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Table 2. Observational information of the XMM- Newton 
data. Column (1): cluster name; Col. (2): Filter in 
wfiicfi cluster was observed; Col. (3): total exposure time; 



Col. (4,5,6) 


: exposure 


time after flare cleaning in 


each de- 


tector. 












Cluster 


Filter t 


exptot (ks) 




^exp (ks) 










MOSl 


M0S2 


PN 


A496 


Thinl 


30.87 


14.44 


14.54 


8.40 


A1689 


Thinl 


39.81 


34.24 


34.64 


27.03 


A2050 


Thinl 


28.41 


24.01 


24.64 


19.87 


A2631 


Thinl 


14.12 


9.32 


11.35 


5.00 


A2667 Medium 


31.16 


22.03 


22.17 


13.88 



majority of the X-ray luminosity. Thus, it is important to 
adequately model the inner parts of a cluster. 

We used the Sersic model and the /5-model to describe 
the surface brightness profiles of our sample with the fits 
presented in Fig|2l For /3-model fits the inner points were 
also included. 

For each cluster an image and the corresponding ex- 
posure maps were created to correct bad pixels and tele- 
scope vignetting. The [0.7-7] keV band is used to derive 
the surface brightness profile, which was obtained adopt- 
ing logarithmic steps between consecutive annuli. Obvious 
point sources or detector gaps were masked and not taken 
into account. We used the STSDAS/IRAF ellipse task. 



3. Gas density distribution 

3.1. Surface brightness profile 

The X-ray emissivity of the intra-cluster medium depends 
strongly on the gas density (ex n^) and weakly on its tem- 
perature (ex T^/^), which means that the observed pro- 
jected emissivity can be used to accurately measure the 
gas density profile. 

In order to compare theory with observations, a 
description of the gas distribution is needed. Using a 
parameter-dependent model for the gas density profile, it 
is possible to reconstruct the 3D X-ray emission of the 
cluster. 

The surface brightness pr ofile has been described for 
many years by the /3-model ("Cavaliere fc Fusco-Femianol 
19761 ). Satellites with better spatial resolution showed a 



significant difference between dat a points and 0-mode\ 
at low radii for cool-core clusters (* J ones fc Formanlll984l : 
Xue fc Wu 2000). An enhancement of the gas density in 
the central region leads to cooling flows and, as a con- 
sequence, the gas density radial profile is cuspy in the 
center instead of having a flat core. Following this obser- 
vational evidence, alternative parametrizations were tried 
to model this em ission excess, notably with the u se of the 
double /3-model (| lke be et al.lll999l : IVikhlinin et~a l. 2006; 
Pratt fc ArnaudI 2005h . In this model, a second term was 
introduced to take into account the change of slope in the 
density profile in the cluster core. This leads to a six pa- 
rameter function {Scooh ^cooU Pcooh Samb, Tamh, Pamb)- 

Another approach towards the descriptio n of the sur- 
face b rightness profile was carried out by Pislar et"al] 
(|l997n . who used the Sersic model to describe the X-ray 
surface brightness profile. 



In a more recent work. lDemarco et al.l (|2003f ) presented 
a detailed analysis of 24 galaxy clusters, showing that the 
Sersic model was adequate to describe the X-ray surface 
brightness profile. Durret et a l. (2005) used the Sersic de- 
scription for Abell 85 and compared it to the /3-model: 
with all points included the Sersic law provided the best 
fit. Moreover, without the inner points, the Sersic and f3- 
model become indistinguishable. 

The central part of cluster is the region where most 
energy is radiated and therefore it is responsible for the 



3.2. Gas density profile modeling 

The surface brightness profile observed is the projection 
of the plasma emissivity and we assume two analytical 
equations to describe the surface brightness profile: the 
/3-model and the Sersic. The former is defined as follows: 



-3/3+0.5 



(1) 



where Tc is the core radius, /3 is the shape parameter and 
So is the central surface brightness. The latter is defined 
by: 



R 



(2) 



where a' is the scale parameter, u (usually represented as 
1/n) is the shape parameter and So is the central surface 
brightness. 

The density profile of the ICM, which emits via ther- 
mal bremsstrahlung, is the deprojection of the surface 
brightness given above and can be given by: 



1+ ( - 



-3/3 



(3) 



for the /3-model and, for the Sersic we have: 



p{r) = po 





-p' 






exp 


i-m 









(4) 



where 
' 2^/^ fourret et al 



a 



1 - 0.6097i/ 

iooi). 



0.05563^^2 ^ 



Note that from the mathematical point of view, the 
Sersic has fewer parameters than the double /3-model and, 
we can adequately compute important quantities such as 
the total mass with no more than three parameters. 

For cool-core clusters, the temperature profile de- 
creases towards the center and, if one assumes that in this 
inner region the derivative of the temperature as a func- 
tion of the radius is positive and that the surface bright- 
ness profile is described by the /3-model we would obtain 
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Fig. 2. Surface brightness profile (black points) of our clusters A496, A2667, A2050, A2631 and A1689. /3-model fits 
are represented by red continuous lines and the Sersic fits are represented by blue dashed lines. The vertical dashed 
lines represent rsoo- For A496 the rsoo is beyond the range of the x-axis. 



a negative mass in the central part. Mathematically, we 
have: 



M(r) (X — r 



d In T d In p 

\ 

d Inr d Inr 



(5) 



which gives the dynamical mass of a cluster. Assuming the 
/3-model, we have: 



d Inp 
lim -— = 0. 

r^o d mr 

For r <C Tc, T oc and: 

d InT 
lim — — = a. 

r^o d mr 

These lead to: 

kTra 



M{r) 



<0, 



G rup II 

which is not a physical solution. 



(6) 



(7) 



(8) 



enough to determine the temperature profile in small bins. 
We excluded from our analysis the most luminous point 
sources which were located inside the FOV of the cameras. 
Each annulus had at least 1200 net counts (i.e., back- 
ground subtracted), which corresponds to a S/N of at 
least 15. 



For the sp ectral fit, we used XSPEC version 11.0.1 
( ArnaudI 1996h and modeled the obtained spectra with 
a Mekal single temperature plasma em ission model 
(jKaastra Mewelll993l : iLiedahl et al.lll995[ ). The free pa- 
rameters are the X-ray temperature {kT) and the metal 
abundance (metallicity). Spectral fits were done in the 
energy interval of [0.7-8.0] keV and, with the hydrogen 
column density fixed at the galactic value (see Table [T]) 
estimated with the tas k nH of FTOOLS (based on 
iDickev fc Lockmanl[T99ol ) 



4. Temperature distribution 

4.1. Temperature profile: spectral analysis 

We extract cluster spectra in concentric annuli around 
the cluster center. On one hand, each annuli had to have 
enough source counts to have good signal-to-noise ratio 
{S/N). On the other hand, the regions must be small 



4.2. Temperature profile modeling and deprojection 



We have chosen an empirical function that can describe 
the overall temperature profiles of cool-core and non cool- 
core clusters: a decrease in temperature towards the clus- 
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Table 3. Fits results. Column (1): cluster name; Col. (2): slope parameter of the /3-model; Col. (3,4): scale radius in 
arcsec and in h^Q kpc of the /3-model; Col. (5): electron number density of the central region calculated by means of 
the /?-model; Col. (6): slope parameter of the Sersic fit; Col. (7,8) scale radius in arcsec and in h^o kpc of the Sersic 
model; Col. 9 electron number density of the central region calculated by means of the Sersic model. 





/3-model 


Sersic 


Cluster 


P 






no 


V 


b1 


a^ 


no 






(arcsec) 


(hro kpc) 


cm~^ 




(arcsec) 


{h~^ kpc) 


cm~^ 


A496 


0.410 ± 0.015 


23.48 ± 0.28 


15.48 =b 0.18 


0.084 


0.390 ± 0.046 


3.6 zb 0.23 


2.37 ± 0.15 


0.036 


A1689 


0.550 ± 0.025 


25.72 ± 0.23 


79.16 ± 0.71 


0.065 


0.650 ± 0.030 


12.78 zb 0.16 


39.32 ± 0.48 


0.012 


A2050 


0.43 ± 0.028 


10.22 lb 0.73 


22.0 ± 1.5 


0.026 


1.285 ± 0.047 


15.33 zb 0.32 


56.4 ± 1.2 


0.010 


A2631 


0.732 ± 0.073 


21.2 ± 2.5± 


82.6 ± 9.6 


0.051 


1.001 ± 0.066 


14.9 zb 1.2 


58.0 ± 4.7 


0.013 


A2667 


0.660 ± 0.017 


30.8 ± 1.9 


113.25 di 6.99 


0.039 


0.570 ± 0.034 


6.3 ± 1.5 


23.1 ± 5.4 


0.022 



ter center, a maximum value and a decrease or a fiat profile 
towards the outskirts. The function is given by Eq. O 









To 


a \ 


h- + l 


V n 


n 



(9) 



^1+1 



where rt is the scale parameter, Tq is the central temper- 
ature and a is the shape parameter. 

The difference between the deprojected temperature 
profile (3D) and the projected (observed) temperature 
profile (2D) fitted above should be less than 10%, 
withi n the observational uncert ainties ( iMarkevitch et al] 
I1999I : iKomatsu fc Seliakl l2QQlh . To be more realistic, 
we deprojected the te mperature profile as described in 



Lima Neto et al.l (|2006f ) to obtain the 3D temperature pro- 



file given by Eq. [TOl 

'd[^'{R)T2Dm 



TsDir) 



r 

Jr 



dR 



OR 



VR^ 



r 

Jr 



dR 



(10) 



dR J ^/W^ 



which will be used in the determination of the dynami- 
cal masses. For non cool-core clusters we simply used the 
mean temperature to infer the total masses. 

To compare the properties of these clusters it is useful 
to determine a representative temperature as the global 
mean temperature. The mean temperature was deter- 
mined extracting cluster spectra inside an external annu- 
lus which encompasses the majority of the X-ray data. For 
each cluster the value of the mean temperature is given in 
Table m 

In Fig. [3] we present temperature profiles with 2D an- 
alytical fits of Eq. O 3D deprojected temperature profiles 
given by Eq. [10] and mean temperatures for clusters in our 
sample. 

5. Dynamical analysis 

The X-ray morphology can give interesting qualitative 
i nsights into the dynamical status of a given cluster 
(| Jones fc FormanI Il984) . Nevertheless, to quantitatively 



Table 4. Best fit results for the parameters of 2D tem- 
perature profiles (Eq.© for A496, A1689 and A2667. 



Cluster 


To 


rt 


a 




(keV) 


(arcsec) 




A496 


1.658 ± 0.069 


726 ± 80 


3.20 ± 0.24 


A1689 


8.16 zb 0.58 


200 zb 26 


-0.27 dz 0.17 


A2667 


3.21 zb 0.39 


184 zb 44 


1.96 ± 0.53 



characterize the dynamical state, density and temperature 
distributions are required. 

During a merger process, the ICM is compressed and 
consequen tly the local density and temperature increase 



(see, e.g.. IVikhlinin et al.ll200ll ). As X-ray substructures 



depend on surface brightness, temperature maps and X- 
ray substructures are auxiliary tools in the analysi s of th e 
dynamical state of clusters (see, e.g.. [Purret et al. 2005 ). 

Star formation efficiency in galaxies is related to the 
dynamical state of galaxy clusters. We present in Fig. [4] 
temperature maps and respective X-ray substructures 
overlaid with the X-ray intensity contours. 

5.1. Temperature maps 

The temperature maps are made in a grid, where 
in each pixel we have made a spectral fit (assuming 
bremsstrahlung + line emission) to determine the temper- 
ature. We set a minimum count number (1200 net counts) 
for a spectral fit with a MEKAL plasma model. If we do 
not reach this minimum count number in a pixel we grad- 
ually enlarge the area up to a 5 x 5 pixel region (or 128'' x 
128'' region). The best fit temperature value is attributed 
to the central pixel. When there are not enough counts 
in the region, the pixel is ignored and we proceed to the 
next one. We compute the effective area files (ARFs) and 
the response matrices (RMFs) f or each region in th e grid. 
This procedure is presented in Durret et al. 

We see from Fig. |4] that temperature maps present 
a level of detail that is lost in the projected temperature 
profiles. For example, hotter blobs are evident even in cool- 
core clusters, which can be globally relaxed clusters. Since 
shock fronts compress and heat the gas during the shocks. 
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temperature substructures can re veal possible h i stories of 



cluster formation, as discussed in lCovone et al. I (|2006al ). 
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Fig. 3. Temperature profile for A469, A1689, A2667, 
A2050 and A2667. For A496, A1689 and A2667 2D an- 
alytical fits given by Eq. [9] are represented by the red 
solid lines, the blue dot-dashed lines represent the com- 
puted deprojected temperature profile given by Eq.fTOland 
the mean temperatures are represented by the horizontal 
dashed lines. For A1689, A2050 and A2631 the mean tem- 
peratures are represented by the horizontal dashed line. 



5.2. X-ray substructures 

In Fig. [4] we present the images of the X-ray substructures, 
which were obtained from the residual images of X-ray P- 
model surface brightness profile fits. The residual image 
is processed using percolation and a filter set to be equal 
to the 3 times the background standard deviation. Only 
pixels with 3a values above the background are going to 
be taken into account to determine the positive (emission 
excess) substructure maps. 

A496 appears to be a very relaxed cluster from its 
X-ray emission contours. However, its temperature map 
shows small regions of substructures, with cooler gas 
seen near the center. This region coincides with the X- 
ray substruc t ures a nd also with the cold front found by 
Dupke et al.l (|2QQ7h . This is evidence of dynamical fossils 
left after a minor merger process. 

Comparing the positive substructure map with the op- 
tical image in Fig. [H we see that for A2050 and A2631 
there is a spatial coincidence between the substructures 
and the bright galaxies, indicating a possible emission 
peak due to star formation heating the surrounding ICM. 
However, comparing the positive substructure map with 
the temperature map, we see that for A496 and A2667 
the positive substructures coincide with the low tempera- 
ture regions. This may be explained by gas infall due to 
the pressure gradients caused by temperature differences, 
consequently increasing the density leading to rise in the 
emissivity. The case of A1689 does not fit in any of the 2 
previously discussed cases, and its over density may simply 
be due to the past motion of halos that compressed the 
ICM, increasing the emissivity. 

From substructure maps we see that all clusters 
but A2667 present evidence of minor mergers. A2667 
could have experienc ed a major merger as proposed by 
Covone et al.l (|2QQ6ah . Even clusters that have a relatively 
smooth appe arance may have und ergone several merger 
events (Durr et fc Lima Netdl2QQ7l ). as revealed by X-ray 
temperature and substructure maps. The deviation from 
the equilibrium hypothes is is n ot significant for the total 
mass estimates (see, e.g.. iRasia et a l. 2006). 

We see that comparing substructure and temperature 
maps with optical images gives us hints about the effects 
of galaxies on the ICM and the effects of ICM motion 
on temperature. In order to investigate any correlation 
between star formation efficiency and cluster dynamical 
state we would need a larger sample, given the intrinsic 
scatter expected for such quantities. 



6. Stellar baryon content 

We can derive the stellar baryon budget by integrating 
the luminosity function and using an appropriate stellar 
mass-to-light ratio. 
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The galaxy luminosity function (GLF) is one of the ba- 
sic statistical properties of galaxies in clusters. The GLF 
gives the number of galaxies per bin of magnitude and in 
most case s, the bright galax i es can be fit by a Schechter 
function ( Schechter fc PressI 1976h . In the case of deep 
observations there appear to be two components in the 
GLF: the bright and faint galaxies. This faint end can be 
seen for magnitudes greater than M^/ = 18 and for these 
cases it is impossible to well describe the GLF with a sin- 
gle S chechter function (iBiviano et al. 19951 : Durret et al 



We have performed analytical fits to the GLF for 
four (A496, A1689, A2050 and A2631) of our five clus- 
ters in order to estimate the total luminosity and then 
the stellar baryon content. For A2667 we use d Lk = 
(1.2 ± 0.1) X lQ^^/i:rn^L ^ given bvlCovone et all (|2QQ6bh . 
based on ISAAC data (jWolfram et al.lll999h . As this lu- 
minosity was computed within R = llO/ifo^kpc in the K 
band, a filter transformation to the Sloan band was done 
(as described in Sec. l6.3p and the total luminosity function 
was extrapolated to rsoo as described in Sect. 16.31 

We calculated the galaxy luminosity function of each 
cluster by means of statistical subtraction. In order to re- 
duce the errors, instead of using all the galaxies we used 
colors to exclude obvious background objects. In Fig. [5] we 
show a {g' — i') x i' color magnitude diagram of the clus- 
ter A2050 showing all galaxies within rsoo- The prominent 
near-horizontal "sequence" is composed predominantly of 
elliptical and lenticular cluster members (cluster red se- 
quence, CRS). For each clusters, we linearly fit the CRS 
{{g' — i') = A^Bx i') and we discarded all galaxies redder 
than {g' — i') = A + 0.3 mag, represented by the solid line 
in Fig. [5l since those are likely background objects. 

The limit in magnitude for each cluster was essentially 
the magnitude where the counts turn over. All galaxies 
that satisfied these two criteria were considered to belong 
to the cluster. 

6.1. Measure of background counts 

In order to accurately compute the cluster luminosity 
function we needed to statistically subtract the back- 
ground from galaxy counts in the cluster direction. Field 
counts should be computed in a region far enough from the 
cluster for it not to be contaminated by cluster galaxies. 
In Fig. [6] we show how the number of galaxies per annulus 
area decreases with radii. From R = 7 rsoo we see that 
galaxy counts becomes flat, meaning that the cluster con- 
tribution is no longer significant and that these counts are 
only due to the background. For this purpose, we selected 
an external annulus whose internal and external radii are 
8 rsoo and 9 rsoo, where we measure the field counts. 
For the number count of galaxies, we have: 



7Vci(^) = A^ci+bkg(m) - 7A^bkg(m), 



(11) 



where Nci{m) is the number of cluster galaxies at a certain 
magnitude m, A/'ci+bkg(^) is the total number of galaxies 
in the line of sight at a certain magnitude, A/'bkg(^) is the 
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Fig. 5. Color- magnitude diagram of galaxies within rsoo 
in the field of view of A2050. The near- horizontal solid line 
indicates the separation between cluster and background 
galaxies. The equation of the solid line is g' — i' = 2.37 — 
0.055i'. The vertical dashed line indicates the i' magnitude 
limit. 
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Fig. 6. Number of galaxies per area of each annulus as a 
function of rsoo for A1689. The background can be deter- 
mined beyond R = Trsoo- 



number of background galaxies at the same magnitude 
and 7 is the ratio between cluster and control field areas. 
After subtracting the background, cluster counts per bin 
of magnitude in the i band were determined. 



6.2. The luminosity functions 

Fitting cluster counts per bin of magnitude with a 
Schechter function described by: 



{m)dm = ^* ioO-4(^*-^)(«+i) exp [-10° 



4(m*— m)i 



(12) 



we obtained the values listed in Tab. [3 where m* is the 
characteristic magnitude and a is the slope of the LF at 
faint magnitudes. 

Schechter function fits were performed using an IDL 
code based on the curvefit function. 
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Table 5. Luminosity function fits results. For A496 we present the results for single and double Schechter fits. 





Schechter 


Cluster 


a 


M* 




A496 


-0.968 ± 0.038 
1.580 ± 0.055 


-22.185 ± 0.099 
-20.551 ± 0.035 


-3.908 ± 0.076 


A1689 


-0.76 =b 0.10 


-21.93 =b 0.19 




A2050 


-1.328 ± 0.072 


-22.22 ± 0.31 




A2631 


-1.19 ± 0.23 


-21.86 ± 0.38 





To obtain m* in absolute magnitude we transform the 
apparent magnitude to absolute magnitude according to: 



where j is the luminosity density of galaxy distribution. 
Assuming a King profile ( King 19781 ) given by: 



m 



2b -b log(^iL/lMpc) + e(z) - K(z), (13) 



where (11 is the luminosity distance, e{z) is the correction 
for evolution and K(z) is the K-correction. Assuming that 
the early-type galaxies are the main population for our 
galaxy clusters at the cluster redsh ift we used K-cor rection 
values for early- types supplied by Poggiantil ( 19971 ). 

For all clusters except A496 we could comfortably fit 
the luminosity functions using a single Schechter function. 
For A496 however, we had to add another power-law com- 
ponent: 



\ -|^Q0.4(m*-m)(Q!2 + l) 



4(m* — m)i 



(14) 



where a and 0^2 are the shape parameters of the bright 
and the faint part of the luminosity function, respectively. 

We did this because the deeper CFHT data for this 
cluster allowed us to sample the faint-end of the luminosity 
function [M'^ > —18) which would not be properly fit 
otherwise. 

6.3. Total luminosity estimation for A2667 

In order to calculate the Sloan i band total luminosity we 
used a filter transformation from K to z band. Assuming 
that the majority of galaxies are ellipticals we convolved 
the spectrum of a z = early-type galaxy taken from 
Coleman etliD (l980) shifted to z = 0.23, with both the 
K and filters, obtaining: 



(17) 



where Tc is the core radius, for which we assumed Tc = 
170 (jGirardi et al.lll995[ ) and jo is the normalization. 
Using the ratio between the total luminosities, L {R = 
r^oo)/L {R = 110 h^o), jo vanishes, and L {R = r^oo) is 
computed as a function of L {R = 110 hj^). 



7. Mass determination 

To estimate the total baryon fraction from the stellar 
components and the ICM, we must choose a fixed ra- 
dius within which the baryonic mass can be computed. 
We choose rsoo, as it is the largest radius for which 
the current X-ray dat a require no model extrapolation 



( Vikhlinin et a 



Lacev fc Cold 



good ) and is about the virial radius 



19931 ). 



Baryonic and total masses are deter mined within this 



radius, which can be derived according to lLima Neto et al 
(^2003,), assuming the isothermal case: 



4.5 X 10^ (3 (kT) 



1/2 



(18) 



where A is the factor by which the mass density exceeds 
the critical density (in this work, A = 500), /3 and Vc are 
the slope and the characteristic radius (given in kpc) given 
by /3-model fit, (/cT) is the mean temperature given in keV 



and /(z; Vl^ 
zf. 



,fiA) = fiA + f^™ (l+2)3-(n™ + l^A-l) (1- 



{i - i^)AB = 1.489. 



(15) 



Using Eq. [151 the Lk given by ICovone et al. ( 2QQ6bl ) 
and the solar mag:nitude in the appropriate bands 
([Fukugita et al.l Il995l : ICo^ I2QQQ[ ) we have Li = 0.98 x 



lO^^/ifo^L© within 110 /ifo kpc 



The total luminosity within r^oo (see Table [7j) can be 
obtained integrating: 



7.1. Gas and total mass determination from 
XMM-Newton data 

Combining the gas density and the deprojected tempera- 
ture profiles, we can derive the total gravitating mass un- 
der the assumption of hydrostatic equilibrium and spher- 
ical symmetry as shown in Eq. O 
The gas mass is given by: 



L 



j 47rr dr. 



(16) Mgas{<r) 



pg Airr'^ dr' 



(19) 
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Fig. 7. Luminosity function for clusters for which we have 
optical data: A496, A1689, A2050 and A2631, respectively. 
For A496 the double Schechter fit (given by Eq. [T4|) is 
represented by the solid line and the single Schechter fit 
(given by Eq.[T2|) is represented by the dotted line. For all 
other clusters we present only the single Schechter fit. 



where pg is the gas density given by Eq. [3] and Eq. [H To 
estimate the central density (po = norrip) we can use the 
normalization of the XSPEC model (i^ parameter): 



14 



EI 



cm 



(20) 



where Da is the angular distance and EI is the emission 
integral, which is model dependent. 

The gas and dynamical mass determinations with both 
models (Sersic and /3 models) are presented in Table[6]with 
the best values in bold. 

For clusters that have a u value of the Sersic fits to 
the surface brightness profile equal to or greater than 1, 
the total mass profiles increase rapidly with radius while 
the gas mass profiles do not have the same behavior. For 
this reason, the mass determination for A2050 and A2631 
is not reliable with the Sersic model and so we do not 
present them in Table [6l 

In order to verify our results, especially for clusters 
that had the surface brightness profile described by the 
Sersic (A496, A1689 and A2667), we compared our es- 
timates of the gas fraction (/gas) presented in Table [6] 
with literature values. Using a dou ble /3-model tode- 



scribe the surface brightness profile, IZhang et al. 
obtained /gas equals 0.102 ± .06 and . 128 ± 0.073 for 



(|2007h 



A1689 and A2667, respectivelv. iLin et al.l (|2003l ) obtained 
/gas = 0.1347 ± 0.32 for A496 also with the double p- 
model. These results show that despite the difference in 
methods, the estimates of the gas fraction are in perfect 
agreement within rsoo- 

7.2. Stellar mass determination from optical data 

To compute the stellar budget of each cluster we used: 



where (^ 



X/sp+(l-/sp)x ( ^ 
blue V 



red 



,(21) 



, is the stellar mass-to-light ratio for late- 

^blue 



type galaxies, (x^)^^^ is the stellar mass-to-light ratio for 
early-type galaxies and /sp is the late-type fraction. 

It has been shown that colors are ^ood estimators 



of morpholofflcal typ e ( Kauffmann et al.ll2003h . Following 
iMateus et all (|2006[ ). we adopted {u - i) = 2.68 to dis- 
tinguish between early- and late- type galaxies, where the 
early- type galaxies are the redder ones. 

From the luminosity function we determined the total 
luminosity according to 



Johi 



L^{L) dL, 



(22) 



where Vd = ^^(^500) is the cluster volume, jobs is the lumi- 
nosity density, which is the total luminosity inside Vd and 
Liim is the luminosity corresponding to the M^/ = — 16 
magnitude. This leads to 



Ltot = VeiL*T[{a + 2l{^)], 



(23) 



where r[a,x] is the incomplete gamma function. 

In order to account for the contribution of fainter 
cluster members, we extrapolate the GLF by integrating 
Eq. [22] with Liim = 0, which is reduced to 



L^{L) dL = L*r(a + 2). 



(24) 
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Using tfiis extrapolation, the total luminosity values would 
be on average 5% greater. Since this extrapolation does 
not greatly change our prior results, we do not use it. 

The late-type and e arly-type stellar mass-to -li^ht ra- 
tios were estimated from Kauffmann et al. ( 2QQ3[ ) and con- 
verted to the i'-band: 

(M/L,). = 0.74 Mo/L©, 



(25) 



for late-types and 
(M/Li), = 1.70 Mq/Lq, 



(26) 



for early- types. 

For A2667 and A496 we do not have enough color in- 
formation to derive the late-type fraction. In these cases 
we assumed that 80% are early-types and 20% are late- 
types to derive a stellar mass-to-light ratio, as these values 
are the average fraction found for the other three clusters. 

Fitting just a Schechter expression for the A496 lu- 
minosity function gives a stellar mass of 1.56 x 10^^ 
instead of 3.51 x lO^^M© obtained from the fit of Eq. [H 
For this cluster the faint-end contributes 33% of the stellar 
mass. This gives us a clue about the error we are introduc- 
ing in not considering the faint-end for the other clusters. 
Our stellar baryon content is underestimated and these 
faint galaxies can contribute differently from cluster to 
cluster. 

7.3. Intracluster light contribution 

Many au thors have d etected a diffuse stellar com- 
ponent (IZibetti et al. 



20051: ICovone et al 



Krick fc BernsteinI 12007 : I Gonzalez et al 




that has possibly been stripped out from galaxies that 
eventual ly merg:ed to form the brightest cluster galaxy 
(BCG) (M urante et aL 2007|) or from the tidal stripping 
of cluster members (|CvDriano et all l2006h . The general 
consensus is that e very cluster has an intracluster light 
(ICL) component (Gonzalez et al.ll2007l ). 

The ICL represents the minor component of the 
baryonic budget, and its contribution to the total lu- 
minosity has been estimated to acc ount for (6 — 22)% 
( Krick fc Bernsteinll200 7). (10-15)% (|Rudick et allEoOG) 
and 10% (IZibettill2007h . Another approach was used by 
Gonzalez et al. who considered in their analysis the 

contribution of the giant BCG and the ICL as a single 
entity (ICL+BCG) that contributes on average 40% and 
30% of the total stellar light, within rsoo and r200 5 respec- 

tively. 

Covone et al. ( 2006al ) studied the diffuse light of 



A2667, but because of the much smaller field-of-view of 
the WFPC2 data they measured only the cD contribution 
to the total light. 

In order to consider the ICL contribution to the baryon 
budget we assumed, for the five clusters in our sample, a 
contribution of ~ 10% to the total luminosity and we used 
the early-type galaxy mass-to-light ratio to obtain the ICL 
mass contribution. 



8. Discussion 

Combining X-ray and optical analysis we show in Fig. [8] 
that for this sample the stellar-to-gas mass ratio is 
anti-correlated with X-ray teni p eratu r e, as proposed 



bv iDavid et all (Il990h 



— - ■■- -Lin et al.l (|2003l ): iGonzalez et al 

(|2007h . iMohr et all (|l999h found that clusters with tem- 
perature below 5 keV have a mean ICM mass fraction sig- 
nificantly lower than that of the hotter clusters, proposing 
an increase of the stellar mass in cooler clusters to main- 
tain the baryon budget. 



0.30 



I 0.15 




0.00 



Fig. 8. MicL+star/^gas as a function of temperature for 
A496, A1689, A2050, A2631 and A2667 (red points). Our 
best fit (Eq. [28|) is repres ented by the solid line. Black 
points were obtained from i Lin et al.l (p003) and the blue 
dot-dashed line was their best fit (Mgtar/^gas kT^-'^^). 

In Fig. [8] we show the stellar-to-gas mass ratio as a 
function of the mean temperature of the clu ster in ou r 
sample and for the clusters studied bv iLin etlil] ( 20031 ). 
There is a clear decrease of Mgtar/^gas with increasing 
gas temperature (which is indicative of the potential well), 
reflecting that the efficiency of converting gas into stellar 
matter is greater in cooler and smaller systems. 

The star formation efficiency in a galaxy can be defined 
in terms of the total stellar mass to gas mass ratio as stated 
above: 



M.as+MiCLy 



MiCL + Msta 



(27) 



Although the number of clusters in this analysis is 
small, a Spearman test of correlation yields p = —0.7 with 
a confidence level of 81%, with all data points included. 
Thus we propose a way to estimate the stellar content 
(diffuse stellar component plus the stars within galaxies) 
of a cluster based on its gas mass and temperature: 



0.18-0.012 



_^gas_ \ 

IkeV / ' 



(28) 
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Table 6. Gas and total mass determination from X-ray data. Column (1): cluster name; Col. (2,3): gas and total 
mass calculated by means of the /^-model; Col. (4,5): gas and total mass calculated by means of the Sersic model; 
Col. (6): gas fraction (Mgas/Mtot) calculated by means of the best model. Values in bold represent the best model for 
the cluster. 





/3-model 


Sersic 




Cluster 


Mgas (lO^'Mo) 


Mtot (lO^'Mo) 




Mtot(lO^"M0) 


/gas 




(< rsoo) 


(< rsoo) 


(< rsoo) 


( < rsoo) 


( < rsoo) 


A496 


1.074 zb 0.040 


1.997 ± 0.012 


0.330 ± 0.040 


3.270 ± 0.19 


0.10 zb 0.014 


A1689 


3.79 ± 0.17 


10.87 it 0.20 


1.672 ± 0.080 


11.14 ± 0.46 


0.150 zb 0.010 


A2050 


0.568 ± 0.079 


4.974 zb 0.094 






0.114 zb 0.016 


A2631 


0.84 ± 0.16 


11.78 ± 0.58 






0.072 zb 0.014 


A2667 


1.789 ± 0.047 


7.06 ± 0.14 


1.132 ± 0.066 


11.63 ± 0.24 


0.0973 zb 0.0060 



Table 7. Results derived from optical data. Column (1): cluster name; Col. (2): total luminosity derived from stars 
within galaxies; Col. (3): fraction of star forming galaxies; Col. (4): stellar mass; Col. (5): ICL mass; Col. (6): stellar 
and ICL mass to gas mass fraction. 



Cluster 


Lstar (1O"'"^L0) 


/. 


Mstar (lO'^Mo) 


MiCL 


Mstar+ICL /Mgas 




(< rsoo) 


(< rsoo) 


(< rsoo) 


(< rsoo) 


(< rsoo) 


A496 


0.0233 ± 0.0021 


0.20* 


0.0351 zb 0.0066 


0.0040 


0.118 lb 0.020 


A1689 


0.0599 ± 0.0023 


0.15 


0.0903 zb 0.020 


0.010 


0.060 zb 0.012 


A2050 


0.0367 zb 0.0010 


0.20 


0.0553 zb 0.0015 


0.062 


0.108 zb 0.040 


A2631 


0.0509 ± 0.0041 


0.20 


0.0768 zb 0.0062 


0.0087 


0.102 ± 0.022 


A2667 


0.160 zb 0.046 


0.20* 


0.240 ± 0.073 


0.027 


0.237 ± 0.029 



* Not measured 



Differently from iLin et al. I Jinni . we found a linear 
correlation between stellar-to-gas mass and temperature, 
which seems to well describe both samples. 

The analysis of the stellar content of A2667 was not 
homogeneous with respect to the other four clusters, re- 
quiring several reasonable assumptions. Its star formation 
efficiency appears in good agreement with the relation de- 
fined by the remaining clusters at the 2cf confidence level. 

We suggest one possible explanation for this trend 
in the context of the hierarchical scenario: during group 
mergers that form rich clusters, different gas galaxy strip- 
ping efficiencies are experienced. In more massive clusters, 
galaxies exhibit higher velocity dispersion and as a con- 
sequence gas stripping is stronger, diminishing star for- 
mation efficiency in galaxies. Thus, physical mechanisms 
depending on virial mass, such as ram-pressure stripping, 
are driving galaxy evolution within clusters. 

Only recently, due to advances in computing capa- 
bilities, cluster simulations have reached adequate lev- 
els of detail to trace star formation. Current numeri- 
cal simulations use a quantity for the cosmic star for- 
mation rate that makes it possible to compare the 
mean efficiency of star formation between halos of dif- 
ferent masses. Hydro dynani ic al simulations performed by 



ior is in conflict with observations, the increa se of the ^as 



fraction with temperature was also found bv iHicks et al 



Springel fc HernquistI (|2003h ; I Romeo et al.l ([2005.) appear 



to confirm that the stellar-to-gas mas s ratio decreas e s with 
increasing temperature. However for Romeo et al.l ([2005') 
this trend is driven by the fact that larger clusters better 
retain the ICM gas, so that the ICM to total mass ratio in- 
creases, while the stellar mass fraction remains fairly con- 
stant with cluster temperature. Although the latter behav- 



(|2007>) who claim that a possible physical explanation for 
lower gas fraction in cooler systems is the higher amount 
of baryonic matter in the form of stars able to mantain 
the baryon budget constant. 

The idea that the actual measurements of /gas in clus- 
ters should be universal is a starting point for inter- 
esting cosmological tests. There has been some debate 
about the dependence of the gas fraction on temperature 
(lArnaud fc EvrardlllQQol : iMohr et al.lll999l : lR.oussel et aT 



2000[ ) and this discussion is still open. Our study suggests 
that the gas fraction is independent of the temperature , 
as shown in Fig. [9] and pointed out in Allen et"al] (|2004l ). 
However, due to the limited sample we cannot draw a firm 
conclusion. To examine in a more systematic and uniform 
way the gas and the stellar fraction variation with tem- 
perature, a larger sample and stellar mass estimations are 
required. 



9. Conclusions 

In this work we combined an X-ray with optical analy- 
sis to study the baryon content (stellar, intracluster stars 
and gas component) of five Abell clusters. Our results are 
robust because of homogenous data processing, a careful 
treatment of the stellar mass-to-light ratio of early and 
late-type galaxies and the use of XMM-Newton data to 
estimate the gas mass. We computed the stellar and the 
gas masses within the same radius. 
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Fig. 9. Gas and total stellar fraction as a function of clus- 
ter temperature. The black points represent the gas frac- 
tion while the red ones are the stellar fraction. For the 
stellar fraction we used both stars in galaxies and the ICL 
contribution. 



From X-rays we derived the surface brightness and 
temperature profiles and we also computed the map of 
substructures for each cluster that gives hints about the 
dynamical state of each cluster. We show that the gas 
fraction obtained by means of the Sersic equation is in 
agreement with values in the literature obtained with the 
double /3-model. 

From optical data we derived the total luminosity and 
the late- type fraction to calculate stellar masses. We as- 
sumed a reasonable contribution of the ICL to the total 
luminosity. 

We confirmed here, usin g a d if jFerent technique , 



early findings b y [ Pavid et all (|l990h : iLin et all ^ 



Gonzalez et al.l ()2QQ7[ ) regarding the depende nce of the 



stellar-to-gas mass ratio on the temperature. iLin et al. 
lave quantified this ratio which agrees with our 

results. 

We showed here that, for this sample of clusters, the 
efficiency of galaxy formation can vary from 6% to 14% 
proving that the star formation efficiency depends on the 
environment, as hydrodynam ical simulation performed by 
Springel fc HernquistI (|2QQ3h have confirmed. 
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